Sufficient conditions for thermal rectification in hybrid quantum structures 
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We analytically identify sufficient conditions for manifesting thermal rectification in two-terminal 
hybrid structures within the quantum master equation formalism. We recognize two classes of recti- 
fiers. In type A rectifiers the contacts are dissimilar. In type B rectifiers the contacts are equivalent, 
but the system and baths have different particle statistics, and the system is (parametrically) asym- 
metrically coupled to the baths. Our study applies to various hybrid junctions including metals, 
dielectrics, and spins. 
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Understanding heat transfer in hybrid structures is 
of fundamental and practical importance for control- 
ling transport at the nanoscale, and for realizing func- 
tional devices Among the systems that fall into 
this category are metal-molecule-metal junctions, the ba- 
sic component of molecular electronic devices 0, and 
dielectric-molecule-dielectric systems, where vibrational 
energy flow activates reactivity and controls dynamics 
Q. Phononic junctions are also captivating and essen- 
tial for understanding the validity of the Fourier's law of 
thermal conduction at the nanoscale 0, 0] • Single- mode 
radiative heat conduction between ohmic metals was re- 
cently detected, showing that photonic thermal conduc- 
tance is quantized [f|. Other hybrid systems with in- 
teresting thermal properties are electronic spin-nuclear 
spin interfaces 0], metal- molecule contacts with exciton 
to phonon energy transfer Q , and metal-superconductor 
junctions @. 

Thermal rectification, namely an asymmetry of the 
heat current for forward and reversed temperature gra- 
dients, has rece ntly attracted considerable theoretical 
0, EL El El EI EE (3 and experimental EL El 
attention. Most theoretical studies, confined to a specific 
realization, have analyzed this phenomenon in phononic 
systems using classical molecular dynamics simulations. 
In this letter we attempt a first step towards an ana- 
lytical understanding of this effect. We establish suffi- 
cient conditions for manifesting thermal rectification in 
a prototype-hybrid quantum model, including a central 
quantized unit (subsystem) and two bulk objects (reser- 
voirs). We identify two classes of thermal rectifiers: (i) 
Type A rectifiers where the terminals are dissimilar i.e. of 
different mean energy (or heat capacity) . (ii) Type B rec- 
tifiers, where the contacts are equivalent, but the reser- 
voirs and subsystem have different statistics, combined 
with unequal coupling strengths at the two ends. We 
manifest that these rectifiers could be realized in several 
subsystems (harmonic and anharmonic) and reservoirs 
(spin, metal, dielectrics), see Fig. [1] 

Consider a 1-dimensional hybrid structure where a cen- 
tral unit Hs interacts with two reservoirs H® {y = L.R) 




FIG. 1: Examples of two hybrid systems treated in this work. 

(a) Single mode heat transfer between a solid and a spin bath. 

(b) Phonon to exciton energy exchange. The central unit can 
represent either a vibrational or a radiation mode. 



of temperatures T v = f3 v 1 via the coupling terms V v , 

H = H L +H° R + Hs + V L + V R . (1) 

The heat current from the left bath into the subsystem is 
given by J L = |Tr ([H° L - H s , V L ]p); (h = 1) [2, where 
p is the total density matrix, and we trace over the system 
and reservoirs degrees of freedom. In steady-state the 
expectation value of the interaction is zero, Tr (^^p) = 
0, and we obtain 



J = iTr[Jp\; J=-[V L ,H S ] 



:[H S ,V R ] 



(2) 



This expression was derived based on the equality Jl = 
—Jr in steady-state. We discard the subscript L in (fj|), 
thus the heat current is defined positive when flowing left 
to right. The system Hamiltonian assumes a diagonal 
form, and we also consider separable couplings 

H s = ^E n \n){n\; 

n 

V y = SB,; S = J2 S m,n\m){n\. (3) 

n,m 

Here S is a subsystem operator and B v is an operator in 
terms of the v bath degrees of freedom. For simplicity 
we set S„ ltn = S ntm . In what follows we consider cases 
where Bl and Br have equal structure but different pref- 
actors. We refer to this case as "parametric asymmetry", 
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or " unequal coupling strength" , rather than " functional 
asymmetry", resulting from dissimilar B's. Note that if 
the commutator [Hs, S] = 0, the heat current trivially 
vanishes. 

We begin and discuss the type A rectifier, constructed 
by adopting reservoirs with distinct properties, as we ex- 
plain below. First we derive a general expression for the 
heat flux in hybrid structures. Formally, if we write the 
heat current as J(T a ,A) = ^ k oik{T a )A k , rectification 
takes place if even terms survive, ak=2n 7^ 0; n = 1, 2..., 
T a = T L + T R and A = T L - T R . 

The initial density matrix p is assumed to be a tensor 
product of system ps and bath ps = Pl{Tl) ® Pr(T r ) 

_H°_ H° 

factors, where p v (T) — e t /Tr„[e t J (y = L,R). 
For convenience we delete the direct reference to time. 
In terms of T„ and A we can write 



1 2T„(ff° +H^)-2A(ff° -H%) 

p(T a ,A) = -e ^ 



1 PS, 



(4) 



where Z = Tt[psPb] is the partition function with the 
trace performed over bath and system degrees of free- 
dom. Furthermore, without loss of generality, the system 
density matrix is assumed to depend initially only on the 
average temperature, ps = Pg (T a /2). 

The expectation value of the energy current is given 
by evaluating J(T Q ,A) = Tr[Jp(T a , A)], where J is an 



Heisenberg representation operator J 
Expanded in powers of A we get 



mt 



J(0)e 



-iHt 



J(T a ,A) 



2A 

J^2 



ao(T a ) 



(2A) S 



ai(T a ) + 0(T a ,A 3 ),(5) 



with the coefficients 
a (T a ) 

ai(T a ) 



Tr 

i T r 
2 



J(H° L -H R )p T \, (6) 

[j(H° L - Hlfp T 
ao(T a )TT[(H° L -H R )p T }. (7) 



Here p T = ^tPl(^-)pr(^-)p s (if) is the density ma- 
trix at the average temperature Tl ^ r with the partition 
function Zt^- Note that the operators are time depen- 
dent, given in their Heisenberg representation. Eq. ([7]) 

was derived using the fact that Tr J px and all its T 
derivatives are zero. Using the definition of the current 
operator ^ we obtain an explicit expression for ao, 



a {T a ) = iTr{[S,H s ](B L Hl + B R H° R )p T }. 



(8) 



Since [Hs, S) ^ [see discussion after Eq. (J3J], the linear 
term in J is finite [201 ] . 

We now examine the onset of thermal rectification, 
i.e. discuss the sufficient conditions for having a± =/= 0. 
While there might be some special values of T a where 
the two terms in ([7J cancel, in general since B v and H® 
are independent operators, the result is finite. Thus, 



in order to manifest rectification it is enough to ana- 
lyze when one of the terms in ([7]) is nonzero, (i) The 
first expression is finite if Tr[Fi?°p T ] ^ Tr[FH R p T }; 
F = [S 7 H s }(B L Hl + B R H R ). The second term is 
nonzero if Tr[iJ^p^] 7^ Ty[H r pt]- Based on (ii) we con- 
clude that rectification emerges when the reservoirs have 
different mean energy, 

(H° L ) = Tr L [p L (T)H° L ] + Tr R [p R (T)H° R ] = (H° R ). (9) 



As an example, consider a bath of 1-dimensional oscilla- 
tors, H® = H^ m + H v ^ 1 , where the kinetic energy H^ m 
is quadratic in momentum and the potential energy per 
particle is C n q n ; n > 2. In the classical limit using the 
equipartition relation we obtain (£f£) = Tl{\+— ). Ther- 
mal rectification thus emerges if the reservoirs have a 
non-identical power n. Note that the separation to three 
segments (L, subsystem, R) is often artificial, as the sys- 
tem can be practically made of a single structure with 
a varying potential energy, e.g. an asymmetrically mass 
loaded nanotube [l7|. Our results manifest that such 
an inhomogeneous structure should rectify heat. Finally, 
we comment that our discussion could be generalized to 
cases where V v = S V B V ; Sl 7^ S R , see Eq. ©. 

We turn to the type B rectifier, and show that for 
equivalent reservoirs rectification emerges when the sub- 
system and reservoirs have different statistics, in conjunc- 
tion with some parametric asymmetry. It is easy to show 
that under ([3]) the steady-state current ([2]) becomes J = 

* J2n,m ^m,nS m ,n^B(BLPm,n), where E m ^ n = E m — E n , 

and Ttb denotes the trace over the L and R reservoirs 
degrees of freedom. Employing the Liouville equation in 
the interaction picture, the elements of the total density 
matrix satisfy 



dt 



= / dr[V(t),[V(r),p(r)]} r 



(10) 



where V = Vl + V R , and V(t) are interaction picture 
operators. Following the standard weak coupling scheme 
21], going to the markovian limit, the heat current re- 
duces to 



J = \zZ E - 



\Sm,n\" Pn x i^n- 



(11) 



where the transition rates are given by k v n ^ m — 
J^ dre l - E "' mT (B v (t)B u (Q)), and the population P n = 
Tre (p n ,n) satisfies the differential equation 



Pn ^ ^ |5m,n| Pjn^fn — >n Pn ^ ^ \Sn 



(12) 



In steady-state P n — 0, and we normalize the popula- 
tion to unity ^ P n — 1. Our description to this point 
is general, as we have not yet specified neither the sub- 
system nor the interfaces. We consider next two rep- 
resentative models for the system Hamiltonian and its 
interaction with the baths. In the first model the sub- 
system is a harmonic oscillator (HO) of frequency u>, 
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Hs = Yin n u\n)(n\. This can describe either a local ra- 
diation mode [6j, 22] or a vibrational mode of the trapped 
molecule (23J. We also take S = Yl n V™\ n )( n ~ M + c - c ) 
motivated by the bilinear form V v oc xi?,,, a; is a subsys- 
tem coordinate [23J. This implies that only transitions 
between nearest states are allowed, 



„v — uv 



dre luJT (B„(t)B v (0)) ; 



ft n-l- 



(13) 



Solving (TT2"]) in steady-state using the rates ([TB"]) , the heat 
current (fTTj) can be analytically calculated, 



;)/fc L 



j)/k R '' 



(14) 



where n v B {uj) — [e' 3 "" — l] 1 is the Bose-Einstein distri- 
bution function at T v = l//3„. 

Our second subsystem is a two-level system (TLS). 
Here Hs = %cr z , and we employ a nondiagonal interac- 
tion S — a x . These terms can represent an electronic spin 
rotated by the environment (jflj. It can also describe an 
anharmonic (truncated) molecular vibration dominating 
heat flow through the junction [l2|, [Hj]. Re-calculating 
the long-time population (|12p. the heat flux reduces to 



j(TLS) 



n^(-uj)/k L +n§ {-uj)/k R 



(15) 



with the rates (|13[) and the spin occupation factor 

rig(u) = [e 13 ^ + l] -1 . Expressions (H]) and (15]) show 
that in the weak coupling limit the effect of the environ- 
ment enters only through the relaxation rates k v , evalu- 
ated at the subsystem energy spacing to. 

We now analyze the general structure of the last two 
expressions, and discuss the onset of thermal rectifica- 
tion. It is clear that if k L (T) = k R (T), i.e. the reservoirs 
and system-bath interactions B v are equivalent, thermal 
rectification is absent as J(A) = — J(— A). On the other 
hand, if k L (T) = f(T)k R (T), resulting e.g. from the 
use of dissimilar reservoirs, the system generally rectifies 
heat besides some special points in the parameter space, 
depending on the details of the model. This case reduces 
to the type A rectifier discussed above. 

However, a more careful analysis of Eqs. (fT4"]) and 
(fTB"]) reveals that rectification prevails if k L (T) = ck R (T), 
c =/= 1 is a constant, given that the relaxation rates' tem- 
perature dependence differs from the central unit particle 
statistics. For example, for a spin (TLS) subsystem we 
require that n s {—uo)/k u = g(T u ), a non-constant func- 
tion of temperature. A TLS asymmetrically coupled to 
two harmonic baths thus rectifies heat [T3|. As we show 
next, the temperature dependence of the rates k u reflects 
the reservoirs statistics. We therefore classify type B rec- 
tifiers as junctions where the system and bath differ in 
their statistics, and the equivalent reservoirs are (para- 
metrically) asymmetrically coupled to the system. We 



specify next the contacts, and exemplify the two classes 
of thermal rectifiers in various hybrid structures. 

Spin bath. Assuming the environment includes a 
set of distinguishable noninteracting spin- 1/2 particles 
(p = 1,2, ..,P), the bath Hamiltonian is given by sum- 
ming all separable contributions, H® = Ylp^lp ano - 
B v = J2 P b v ,p- The relaxation rate (fTB"]) reduces to 



1] 

2 

5{oj + 



(16) 



and T u (uj) 
e p (0) - e p (l)), with 



the 



where 

p-particle eigenstates \i) and eigenvalues e p (i) (i = 0, 1). 

Solid/ Radiation field (harmonic bath). This bath in- 
cludes a set of independent harmonic oscillators, cre- 
ation operator a u j . System-bath interaction is further 

assumed to be bilinear, 



y., -v. 



",3 



), where X„ j are the system-bath cou- 



pling elements. This leads to the relaxation rate (fTB"]) 



k v 



-Y v B (w)n v B (-w), 



(17) 



where T b (oj) — 2ir J^j \^u,j | 2 — ^) is an effective 
system-bath coupling energy and n v B is the Bose-Einstein 
distribution function at temperature T v . 

Metal. As a final example the contact is made metallic, 
including a set of noninteracting spinless electrons, cre- 
ation operator c v i} The bath operator coupled to the sys- 
tem allows scattering between electronic states within the 
same lead, H° = V, < ,'V.,'V.<: B v = J2i,j v ^i,j c l,i c ",j- 
The transition rate (JTBJ can be written as 



k v = -2im B {-uj) \v u ,i,j\ S(ej - e + u)) 

id 

x[rip(ei) - n F (e l + u)], 



(18) 



with the Fermi-Dirac distribution function n F (e) = 
[ e /Me-/*>0 + 1]- 1 a t the chemical potential fj, u . One could 
also write 



r = -n v B (-u)A"(T V) w) ) 



(19) 



2ir J de [n F (e)-n F {e + u)]F v {e). 



where k v (T v u) 

The function F v {e) = \ v v,i,j | z <5(e — Cj + oj)8(ei — e) 
depends on the system-bath coupling elements and the 
specific band structure. Assuming that the density of 
states slowly varies in the energy window uj, this func- 
tion could be expanded around the chemical potential 
[22l ]. If the Fermi energy is much bigger than the con- 
duction band edge we obtain K v (T v w) sw T F (l + Sv^f-j , 

where T F = 2irLuF u (fi„), and S is a constant of order one, 
measuring the deviation from a flat band structure near 
the chemical potential [22| . 

Consider for example a type A solid-HO-metal recti- 
fier, representing an electronic to vibrational energy con- 
version device, see Fig. Gib). This system may be real- 
ized by attaching an insulating molecule to a metal (STM 
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FIG. 2: Spin-HO-Metal rectifier (full) and a Spin-TLS-Metal 
rectifier (dashed). Rectification ratio is presented as a func- 
tion of the subsystem spacing. T a — 0.5, A = 0.1, Sr — 0.2, 
fiR = 1, Main plot: The subsystem is equally coupled to the 
two ends, Tl = = 1. Inset: The effect can be tuned by 
manipulating system-bath couplings, Tl — 1, Tp = 0.05. 



tip), while the underneath surface is insulating. Setting 
the coupling strength at both contacts to be the same, 
uj) — Tp(oj), we can calculate the rectification ratio 



1 B 

n = 



J(-A) 



using dHJ), (HH and HH) 



n ii R i 



(1 + *h^)(2- 



1-A-^. (20) 



Therefore, if the metal density of state varies with energy 
(5 =/= 0), thermal rectification is presented [22]. Interest- 
ingly, we can show that in a solid-TLS-metal junction 1Z 
could be modulated to be greater or smaller than one by 
varying the gap u>. Thus, phonon-to-exciton heat conver- 
sion can be made effective, while the exciton-to-phonon 
route becomes ineffective, and vice versa. As a second- 
type B- example consider a solid-HO-solid structure. In 
the classical limit [see ([TBI and (ITT)) ] we get 



J(T a ,A) 



2T%T L + 2T«T R 

A / A 
— 1 - x 1- a 

T n \ T„ 



T 2 



(21) 



with x = (r£ - rf )/(r| + rf ). 

that even A terms, i.e. thermal rectification, are directly 
linked to the parametric asymmetry. Fig. [2] further dis- 
plays the tunability of a spin-subsystem-metal junction. 
In the classical limit {uj < T„) rectification can be sub- 
stantial, while in the quantum regime the effect is sup- 
pressed. Modifying the system-metal coupling strength 
largely controls the rectification ratio (inset). 

To summarize, while previous studies were focused on 
a specific realization, typically limited to the classical 
regime, based on numerical simulations, we have analyti- 
cally deduced sufficient condition for the onset of thermal 
rectification in generic hybrid structures at the level of 
the quantum master equation: (i) The reservoirs should 
be made dissimilar, e.g. rectification emerges in an an- 
harmonic junction where the L and R segments have 
different potentials, (ii) The contacts could be of the 
same type, but their statistics should differ from that 
of the system, combined with some parametric asymme- 
try, e.g. a boson-spin-boson junction rectifies heat when 
^ Tg. Our study applies to various interfaces: met- 
als, insulators and noninteracting spins. The central unit 
could represent a radiation mode, a vibrational mode, or 
an electronic excitation. 

Anharmonic interactions were in particular pointed 
out responsible for thermal rectification. While it is obvi- 
ous that perfectly harmonic systems cannot rectify heat 
[H, not all anharmonic-asymmetric systems do bring in 
the effect. Consider for example a three-segment nonlin- 
ear oscillators chain where all units have identical poten- 
tials, but the central part is asymmetrically connected to 
the two terminals. According to our analysis this system 
will not rectify heat. 

Our study manifests control over energy transfer at the 
nanoscale, important for cooling electronic and mechan- 
ical devices and for controlling molecular reactivity. We 
also demonstrated that thermal rectification is an abun- 
dant effect that could be observed in a variety of systems, 
phononic [ItJ , electronic [H[ , and photonic [1, [12] ■ 
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